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A microscopic approach has been employed to study the kaon productions in heavy ion collisions. 
The momentum integrated Boltzmann equation has been used to study the evolution of strangeness 
in the system formed in heavy ion collision at relativistic energies. The kaon productions have been 
calculated for different centre of mass energies (^snn) ranging from AGS to RHIC. The results have 
been compared with available experimental data. We obtain a non-monotonic horn like structure for 
/tv^ when plotted with ^snn with the assumption of an initial partonic phase beyond a certain 
threshold in y^ijvjv. However, a monotonic rise of j-K^ is observed when a hadronic initial state 
is assumed for all -^/snn. Experimental values of K~ /-k" are also reproduced within the ambit of 
the same formalism. Results from scenarios where the strange quarks and hadrons are formed in 
equilibrium and evolves with and without secondary productions have also been presented. 

PACS numbers: 25.75.-q,25.75.Dw,24.85.+p 



I. INTRODUCTION 

The lattice simulation of Quantum Chromody- 
namic equation of state (EoS) predicts that the prop- 
erties of nuclear matter at extreme densities and/or 
temperatures are governed by the partonic degrees 
of freedom Ih-Q- ^ series of experiments have been 
performed [J| and planned [5| to produce such a 
partonic state of matter, called Quark Gluon Plasma 
(QGP) by colliding nuclei at ultra-relativistic ener- 
gies. Rigorous experimental and theoretical efforts 
are on to create and detect such a novel state of 
matter 0]. Various signals have been proposed for 
the detection of QGP - the pros and cons of these 
signals are matter of intense debate. The study of 
the ratio, = ivr+/7r+ is one such currently de- 
bated issue, i?"*" is measured experimentally |7HlO{ 
as a function of centre of mass energy {^/snn)- It 
is observed that the i?"*" increases with ^/sWn and 
then decreases beyond a certain value of ^ysNN giv- 
ing rise to a horn like structure, whereas the ratio, 
= /tt^ increases faster at lower i/J/vw and 
tend to saturate at higher ^/snn- 

Explanation of this structure has ignited intense 
theoretical activities (TT| - [l9| . Several authors have 
attempted to reproduce the K'^ /n^ ratio using dif- 
ferent approaches. While the authors in use a 
hadronic kinetic model, in Ref. [13j high mass un- 
known hadronic resonances have been introduced 
through Hagedorn formula to describe the data. In 
Ref. |14| a transition from a baryon dominated sys- 
tem at low energy to a meson dominated system at 
higher energy has been assumed to reproduce the ra- 
tio /tt^. The release of color degrees of freedom 
is assumed in fll| beyond a threshold in ^snn (re- 
sulting in large pion productions) or the production 



of larger number of pions than kaons from higher 
mass resonance decays has also been employed [l^ 
to explain the data. In the present work we employ 
a microscopic model for the productions and evo- 
lution of strange quarks and hadrons depending on 
the collision energy. Here we examine whether the 
/v +/7r+ experimental data can differentiate between 
the following two initial conditions or two scenarios 
- after the collisions the system is formed in: (I) 
the hadronic phase for all ^snn or (II) the par- 
tonic phase beyond a certain threshold in -^snn- 
Other possibilities like formation of strangeness in 
complete thermal equilibrium and evolution in space 
time (III) without and (IV) with secondary produc- 
tions of quarks and hadrons have been considered. 
(V) Results for an ideal case of zero strangeness in 
the initial state has also been presented. 

In the context of strangeness enhancement as sig- 
nal of QGP formation, similar approach i.e. the as- 
sumption of an initial state where non-strange sec- 
tors are in equilibrium but the strange degrees of 
freedom are out of equilibrium (having density much 
below their equilibrium values) were considered in 
the 1980's. The strangeness production in a decon- 
fined (partonic) phase is enhanced compared to the 
their production in the confined (hadronic) phase 
primarily because even the lightest strange hadrons, 
the kaons are much heavier than the strange quarks. 
Moreover, the strange quark has more degrees of 
freedom (six) in a deconfined matter compared to 
kaons. Therefore, the strangeness production dur- 
ing the space time evolution of the system for par- 
tonic initial state will be enhanced compared to the 
hadronic initial state, hence the enhanced produc- 
tion of strangeness could be an efficient signal for 
deconfinement [20l - [2^ . In contrast to these stud- 
ies Gazdzicki and Gorenstein [111, HI, HJ] within the 
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ambit of statistical model considered the strangeness 
production where both the strange and non-strange 
degrees of freedom are in thermal equilibrium and 
the production of strangeness during the expansion 
stage is ignored. In the present work we would like 
to compare the results on kaon to pion ratio from 
these two contrasting scenarios. 

We assume that the non-strange quarks and 
hadrons are in complete thermal (both kinetic and 
chemical) equilibrium and the strange quarks and 
strange hadrons are away from chemical equilibrium. 
Therefore, the evolution of the strange sector of the 
system is governed by the interactions between the 
equilibrium and non-equilibrium degrees of freedom. 
The momentum integrated Boltzmann equation pro- 
vides a possible framework for such studies. Simi- 
lar approach has been used to study the sequential 
freeze-out of elementary particles in the early uni- 
verse [2^. 

For the strangeness productions in the partonic 
phase we consider the processes of gluon fusion and 
light quarks annihilation. For the production of 
Js'+ and K~ an exhaustive set of reactions involving 
thermal baryons and mesons have been considered. 
The time evolution of the densities are governed by 
the Boltzmann equation. 

The paper is organized as follows. In the next sec- 
tion the rate of strangeness productions in the par- 
tonic and hadronic phases are discussed. The space 
time evolution of the system is presented in section 
III. Results are presented in section IV and finally 
section V is devoted to summary and conclusion. 



II. STRANGENESS PRODUCTIONS 

The productions of s and s in the QGP and the 
and in the hadronic system are discussed 
below. 



A. Strange quark productions in the QGP 
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FIG. 1: Rate of production of s quark from gg — ^ ss and 
qq — )■ ss with temperature. 
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FIG. 2: The rate of kaon production from dominant me- 
son meson interactions with temperature. 



The two main processes for the strange quark pro- 
ductions are gluon fusion {gg — ^ ss) and quark((7)- 
antiquark (g) annihilations {qq ss). The cross 
sections in the lowest order QCD is given by [20|: 

87ra? , 2m? , , , 

Cjqq^,s = ^{l + —)W{S) (1) 



and 



2-1:01 . 1 -1 , \ 7 31r7i^ 
0-93^35 = ^ — [G{s)ianh w{s) - - ^ — w{s)\ 

oS o oS 

(2) 

where m is the mass of strange quark, s = {pi +P2)'^, 
is the square of the centre of mass energy of the col- 
liding particles, pi are the four momenta of incoming 



particles, G ~ 1 + Arri^ / s + ni'^ / s'^ , w(s)=(l — 4to^/s) 
and as is the str ong coupling constant that depends 
on temperature (26|. 



B. if^ and K productions in the hadronic 
system 

The rate of K^{us) and K~{us) productions in 
the hadronic phase can be categorized as due to (a) 
meson- meson (MM), (b) meson-baryon {MB) and 
(c) baryon- baryon {BB) interactions. In the present 
paper we quote only the main results for kaon pro- 
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ductions in the hadronic matter and refer to 27] for 10 

details. jq"' 

(a) For the first category MM KK, we consid- _g 
ered the following channels: rnr -> KK, pp — > KK, 

np — )■ KK* and np — )■ K*K. The invariant ampli-rr i(f^ 

tude for these processes have been calculated from *^ lo 

the following Lagrangians [13] ■ For the K*Ktt ver- ^ ^ 




tex the interaction is given by, 

Ck'K. = gK'K.K*^T[K{d,,n) - )7r] 

Similarly for the pKK vertex the interaction is, 

CpKK = gpKK[KT{d^K) - {d^K)TK]p^' (4) 

The isospin averaged cross section [a) for A/M — S' 
ii'A' (i.e.,_7r7r ^- pp and Trp i^TiC*, 

7rp — )• K* K) is evaluated by using. 
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where P and P' are the three-momenta of the meson 
and kaons in the centre-of-mass frame, x is the cosine 
of the angle between P and P' . M{s,x) is the iso- 
spin averaged squared invariant amplitude. 

(b) For meson baryon interactions the dominant 
channels are: n N — > A if, pN — > AK, tt N — >• 
N K K and tt N — > N tt K K. The isospin averaged 
cross section is given by jisj j: 



f^MB^YK — 
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Ji, Fi and rrii are the spin, width and mass of the 
resonances, (25 -I- 1) is the polarization states of the 
incident particles, k is the centre of mass momentum 
of the initial state. P™ and B°^* are the branching 
ratios of initial and final state channels respectively. 
The index i runs over all the resonance states. For 
interactions ttN — >■ AK, pN — ^ AK we have con- 
sidered 7V*(1650), A^2*(1710) and 7V3*(1720) as the 
intermediate states. Values of various hadronic 
masses and decay widths are taken from particle 
data book 

(c) For the last category of reactions i.e. for 
baryon baryon interactions [29l - l3l| the dominant 
processes are: NN N AK, 7VA_-> N AK, 
A A NAK, NN N N KK, NN 
NNrnrKK and NN^NNnKK. 

The isospin averaged cross section of kaon produc- 
tion from the process like N1N2 — N3AK is given 
by MM 
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FIG. 3: Rate of kaon productions from the meson-baryon 
interactions with temperature. 
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pion is the intermediate particle for the above in- 
teraction, niN is the mass of N, W is the total 
energy in the centre of mass system of pion and 

^2, Wmin = rriK + mA,Wmax = S^^^ " ™7V ■ 

2™^ — 2EE' ± 2pp' where p, p' are the mo- 
menta and E, E' are the energies of A^i and N3 
respectively. We take Jt^nn =1 and to constrain 
the finite size of the interaction vertices we use 
the form factor F — (A^ — mJ)/(A^ — q^)- ctq is 
the isospin averaged cross section of ttN2 — > AK. 
Cross sections for the processes: A^A — >■ NAK and 
A A — NAK have been taken from [sO]. The cross 
section of other reactions e.g. N N N N K K , 
NN ^ N NnnK K ernd N N ^ N N tt K K have 
been taken from (3]| . In a baryon rich medium, 
K~ gets absorbed due to its interaction with the 
baryons. The reactions K~p Avr", K^p — > an^, 
K~n — >■ ap, K~p — )■ K^n, K^n — >■ K^n have 
been considered for K^ absorption (sij in the nu- 
clear matter. 



C. Rate of strangeness productions 

dN/d'^x{= R), the number of s quarks produced 
per unit time per unit volume at temperature T and 
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FIG. 4: Rate of kaon productions from baryon baryon FIG. 6: Comparison between rates of kaon productions 
interactions with temperature. from MM and MM + MB interactions with temperature. 




FIG. 5: Rate of Kaon productions from mcson- 
meson(MM) interactions, meson-baryon(MB) and 
baryon-baryon(BB) interactions at different center of 
mass energies 



baryonic chemical potential /is is given by 

where piS are the momenta of the incoming parti- 
cles and fipiYs are the respective phase space dis- 
tribution functions (through which the dependence 
on T and are introduced), Vrei = |wi — W2I is the 
relative velocity of the incoming particles and a is 
the production cross sections for the reactions. The 
same equation can be used for kaon production by 



FIG. 7: Total if^ and K production rates with temper- 
ature at center of mass energy=7.6 GeV and 200 GeV. 



appropriate replacements of phase space factor and 
cross sections. 



III. EVOLUTION OF STRANGENESS 

The possibility of formation of a fully equilibrated 
system in high energy nuclear collisions is still a 
fiercely debated issue because of the finite size and 
life time of system. In the present work we assume 
that the strange quarks or the strange hadrons (de- 
pending on the vahic of ^snn) produced as a result 
of the collisions are not in chemical equilibrium. The 
time evolution of the strangeness in either QGP or 
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hadronic phase is governed by the momentum in- 
tegrated Bohzmann equation. We have assumed 
that the initial density of strange quarks or kaons 
(depending on the initial conditions (I) or (II)) is 
20% away from the corresponding equilibrium den- 
sity. We will comment on the amount of deviations 
from chemical equilibrium later. 



A. Evolution in QGP and hadronic phase 

The momentum integrated Boltzmann equation 
has been applied to study the freeze-out of elemen- 
tary particles during the thermal expansion of the 
early universe [25|. In the present work we fol- 
low similar procedure to study the evolution of the 
strange quarks and anti-quarks in the QGP phase or 
kaons in the hadronic phase. The coupled equations 
describing the evolution of i (particle) and j (anti 
particle) with proper time r is given by: 
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T 



^. (9) 



where, Ui (rij) and Ui'^'^ i'l^j'^'') ^ire the non- 
equilibrium and equilibrium densities of i (j) type 
of particles respectively. Ri is the rate of produc- 
tion of particle i at temperature T and chemical 
potential /is, r is the proper time. First term on 
the right hand side of Eq. |9] is the production term 
and the second term represents the dilution of the 
system due to expansion. The variation of tempera- 
ture and the baryonic chemical potential with time 
is governed by the hydrodynamic equations (next 
section) . The indices i and j in EqlU] are replaced by 
s, s quark in the QGP phase and by K^,K~ in the 
hadron phase respectively. 



B. Evolution in the mixed phase 

For higher colliding energies i.e., ^/s > 8.76 GeV 
an initial partonic phase is assumed. The hadrons 
are formed at a transition temperature, Tc = 190 
MeV through a first order phase transition from 
QGP to hadrons. The fraction of the QGP phase 
in the mixed phase at a proper time t is given 
by Hi, HI: 



-(r 1) 

r — 1 T 



(10) 



where tq (th) is the time at which the QGP (mixed) 
phase ends, r is the ratio of statistical degeneracy in 



QGP to hadronic phase, 
are governed by |32l |: 



The evolution of the kaons 
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Similar equation exist for the evolution of s and s 
quarks in the mixed phase (see [s^ for details). In 
the above equations ///(r) = 1 — /q(t) represents 
the fraction of hadrons in the mixed phase at time 
T. The last term stands for the hadronization of s{s) 
quarks to K^{K~) t34j. Here (5 is a parameter 
which indicates the fraction of s{s) quarks hadroniz- 
ing to K^{K~). 5 = 0.5 indicates the formation of 
A' + and in the mixed phase because half of the 
s form _ftr+ and the rest hadronize to . 



C. Space time evolution 

The partonic/hadronic system produced in nu- 
clear collisions evolves in space-time. The space- 
time evolution of the bulk matter is governed by the 
relativistic hydrodynamic equation: 







(12) 



with boost invariance along the longitudinal direc- 
tion [111. In the above equation T'"" = {e+P)u^'u''- 
g^'^ P, is the energy momentum tensor for ideal fiuid, 
e is the energy density, P is the pressure and is the 
hydrodynamic four velocity. The net baryon number 
conservation in the system is governed by: 



5,.(nBM^) = 



(13) 



where is the net baryon density. Egs. [T^and[T51 
have been solved (see [H, [13] for details) to ob- 
tain the variation of temperature and baryon den- 
sity with proper time. The initial temperatures cor- 
responding to different ^snn are taken from Table 
I. The baryonic chemical potential at freeze-out are 
taken from the parametrization of /is with ^snat 
[HI (see also [l3|) and the baryonic chemical po- 
tential at the initial state is obtained from the net 
baryon number conservation equation. 

The initial temperatures of the systems formed 
after nuclear collisions have been evaluated from the 
measured hadronic multiplicity, dN/dy by using the 
following relation: 



2^4 
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(14) 
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where C(3) denotes the Riemann zeta function, R 
is the transverse radius l.l{Npart/'2y^^ , Npart is 
the number of participant nucleons ] of the coUiding 
system , ti is the initial time and gef / is the statis- 
tical degeneracy. Initial temperatures for different 
y/sNN are tabulated in table I. 

TABLE I: Initial conditions for the transport calculation. 
Colliding energies are in centre of mass frame 
VsNN (UeV) T, (GeV) Tc (GeV) 



3.32 


0.115 




3.83 


0.128 




4.8 


0.150 




6.27 


0.160 




7.6 


0.187 




8.76 


0.210 


0.190 


12.3 


0.225 


0.190 


17.3 


0.25 


0.190 


62.4 


0.3 


0.190 


130 


0.35 


0.190 


200 


0.40 


0.190 



IV. RESULTS AND DISCUSSION 

The variation of the number of strange anti-quarks 
produced per unit volume per unit time with tem- 
perature has been displayed in Fig. [1] for a bary- 
onic chemical potential = 107 MeV. It is ob- 
served that the process of gluon fusion dominates 
over the qq annihilation for the entire temperature 
range under consideration, primarily because at high 
A*b(= 3/ig) the number of anti-quarks is suppressed. 
In Fig. [21 the production rate of from the 
MM — KK type of reactions has been depicted for 
^snn = 7.6 GeV. The production rate from pion 
annihilation dominates over the reactions that in- 
volves p mesons, because the thermal phase space 
factor of p is small due its larger mass compared 
to pions and smaller production cross section. Re- 
sults for interactions involving mesons and baryons 
are displayed in Fig. [3] It is observed that the in- 
teractions involving pions and nucleons in the initial 
channels dominate over that which has a p meson 
in the incident channel. In fact, contributions from 
the reactions p N ^ A K has negligible effect on 
the total productions from the meson baryon inter- 
actions. The kaon production from baryon-baryon 
interaction is displayed in Fig. 21 The contributions 
from A — > N AK dominates over the contribu- 
tions from N N N AK and A A ^ iV A X for 
the temperature range T — 120 to 180 MeV. 

In FigOthe rates of productions from meson- 
meson interactions has been compared with those in- 
volving baryons i.e. with meson-baryon and baryon- 
baryon interactions for different ^snn (different 



/is). The results clearly indicate the dominant role 
of baryons at lower collision energies which dimin- 
ishes with increasing ^snn. At low temperature the 
baryonic contribution is more than the mesonic one 
for lower beam energy. Rate of productions (from 
MB-I-BB interactions) at ^snn=4.8 GeV is more 
compared to the rates at ^snn=7.6 and 200 GeV, 
since p,B at ^snn=4.8 GeV is more (see table-II). 
Production rate from pure mesonic interactions does 
not depend on pB hence same for all. It is quite 
clear from the results displayed in FigjSl that more 
the baryonic chemical potential (lower the centre of 
mass energy), more is the rate from BE and MB in- 
teractions compared to MM interactions. For a sys- 
tem having lower chemical potential (higher centre 
of mass energy) the rate of production from mesonic 
interactions is dominant. A comparison is made 



TABLE II: Chemical potential for different centre of 
mass energies 



Center of mass energy(-ysjvjv) 
(in A GeV) 


Chem. potential (pb) 
(MeV) 


3.32 


595 


3.83 


568 


4.8 


542 


6.27 


478 


7.6 


432 


8.76 


398 


12.3 


321 


17.3 


253 


62.4 


86 


130 


43 


200 


28 



between rates of kaon productions from meson me- 
son interactions (MM) and meson-baryon (MB) plus 
baryon-baryon (BB) interactions for ^snn = 7.6 
GeV. At this energy baryons and mesons are equally 
important as shown in Fig. [51 

In Fig. [Tithe net rates of productions for and 
K~ have been depicted for ^snn = 7.6 GeV (left 
panel) and 200 GeV (right panel). At -/snn — 7.6 
GeV the production of dominates over K~ for 
the entire temperature range. However, for large 
y/fiNN (low Pb) the productions of and are 
similar. The strong absorption of the by nu- 
cleons in a baryon rich medium resulting in lower 
production yield of K~ compared to K^. This 
may be contrasted with the experimental findings 
of BRAHMS experiment [3^ where it is observed 
that at mid-rapidity (small pB due to nuclear trans- 
parency at RHIC energy) the and K~ yields are 
similar but at large rapidity (large pb) yield is 
smaller than K'^ due to large K" nucleon absorp- 
tion. 

In Fig. [51 the variations of i?+ with ^snn are de- 
picted. The experimental data on is well re- 
produced if a partonic initial phase (scenario-II) is 
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assumed beyond ^snn=8.7 GeV. A "mindless" ex- 
trapolation of hadronic initial state (scenario-I)for 
all the y/sfm up to RHIC energy show an increasing 
trend in disagreement with the experimental data 
at higher ^snn- In both the scenarios, I and II, 
the curves at higher ^snn (RHIC energies) becomes 
flatter. That is because at higher energies the 
productions in the hadronic phase are dominated by 
mesonic interactions and the production rates from 
mesons are same for all ^snn for a given tempera- 
ture range. But at lower energies the rates of kaon 
productions are dominated by the effective interac- 
tions among the baryonic degrees of freedom. The 
composition of matter formed in heavy ion collision 
changes from a matter dominated by baryons to a 
matter dominated by mesons with the increase in 
colliding energy. The hb changes from 86 MeV to 
28 MeV as y^s^ varies from 62.4 GeV to 200 GeV 
(Table II). The change in the production in the 
hadronic phase due to the change in /i^ mentioned 
above is marginal - resulting in the flatness in i?+ at 
higher energies. The decrease of the value of the i?"*" 
beyond ■y/SNN=7.6 GeV showing 'horn' like struc- 
ture happens only when an initial partonic phase 
is considered. Such a non-monotonic behaviour of 
i?+ can be understood as due to larger entropy pro- 
ductions from the release of large colour degrees of 
freedom (resulting in more pions yield) compared to 
strangeness beyond energy 7.6 GeV . 

In Fig. ini the variations of R~ with ^^snn is dis- 
played. R~ has a lower value compared to at 
lower energies since K~ get absorbed in the bary- 
onic medium. At higher energies is closer to 
because production of and _ftr~ is similar 
in baryon free medium, which may be realized at 
higher collision energies. 

In Fig. [To] the i?+ is depicted as a function of 
Y^SNTV for other scenarios (III, IV and V). on the 
one hand when the strange quarks and kaons are 
formed in complete equilibrium but their secondary 
productions are neglected during the evolution (sce- 
nario III) then the data is well reproduced. On the 
other hand, in the scenario (IV) when the system is 
formed in equilibrium (as in III) but the productions 
of strange quarks and kaons are switched on through 
secondary processes then the data is slightly overes- 
timated at high ^SNTV- However, we have seen that 
the data is also reproduced well in the scenario II as 
discussed above. This indicate that the deficiency 
of strangeness below its equilibrium value as consid- 
ered in (II) is compensated by the secondary pro- 
ductions. In scenario V we assume that vanishing 
initial strangeness and observed that the production 
of strangeness throughout the evolution is not suffi- 
cient to reproduce the data. The productions from 
secondary processes are small but not entirely neg- 
ligible (V). In Fig. [11] the R~ has been displayed as 



0.45 

0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 



• KVti^ Data 

- Scenario-I 

— Scenrio-II 




10 

1/2 



100 



(GeV) 



FIG. 8: /■k'^ ratio for different centre of mass en- 
ergies. Scenario - I represents for pure initial hadronic 
scenario for all centre of mass energies. Scenario - II 
represents for the calculation with hadronic initial con- 
ditions for low y/SNN and partonic initial conditions for 
higher y/SNN- See text for details. 
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FIG. 9: K~ /n~ ratio for different centre of mass ener- 
gies. Scenario-I represents for pure initial hadronic sce- 
nario for all centre of mass energies. Scenario-II repre- 
sents for the calculation with hadronic initial conditions 



for low y/sJfN and partonic initial conditions for higher 
y^sjvjv. See text for details. 



a function of ^snw. A trend similar to the results 
shown in Fig. [TU] is observed. The data is overesti- 
mated for the intermediate ^snw in the scenario IV, 
reproduced well in scenario HI and underestimated 
for the scenario V. 
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FIG. 10: K'^ /■k'^ ratio for different centre of mass en- 
ergies. Scenario-Ill assumes complete equilibrium of 
strange quarks and hadrons. The production through 
secondary processes have been ignored. Scenario IV is 
same as III with secondary productions processes are on 
and scenario V represents zero strangeness initially but 
secondary productions are switched on. 
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FIG. 11: Same as Fig.[TO]for K' /tt- 



V. SUMMARY AND CONCLUSIONS 



The evolution of the strangeness in the system 
formed in nuclear collisions at relativistic energies 
have been studied within the framework of momen- 
tum integrated Boltzmann equation. The Boltz- 
mann equation has been used to study the evolution 
of s and s in the partonic phase and K~ and in 
the hadronic phase. The calculation has been done 
for different centre of mass energies ranging from 
AGS to RHIC. We get a non-monotonic variation of 
A'+/7r+ with y/SNN when an initial partonic phase 
is assumed for ^/SNN = 8.76 GeV and beyond. A 
monotonic rise of /tt^ is observed when a pure 
hadronic scenario is assumed for all centre of mass 
energies. The K~ /Tr~ data is unable to differentiate 
between the two initial conditions mentioned before. 



Some comments on the values of the initial param- 
eter are in order at this point. We have seen that 
a 10% variation in the initial temperature does not 
change the results drastically. We have assumed that 
the initial density of strange quarks or kaons depend- 
ing on the scenario (I) or (II) is about 20% away from 
the corresponding equilibrium density. Results from 
a scenario where strange quarks or kaons are formed 
in complete equilibrium and the production is ig- 
nored during the evolution then the data is well re- 
produced (scenario III). If the the strangeness is pro- 
duced in equilibrium and the production is included 
during the expansion stage then the data is overes- 
timated. However, if the system is formed with zero 
strangeness then the theoretical results underesti- 
mate the data substantially. This indicate that the 
production of strangeness during the expansion of 
the system is small but not entirely negligible. The 
deficiency assumed in scenario II is compensated by 
the production during evolution. 
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